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Nonlinear  Gulf  Stream  interaction  with  the  deep 
western  boundary  current  system:  observations 
and  a numerical  simulation 

By  David  E.  Dietrich,  Avichal  Mehra  f,  Robert  L.  Haney  $,  Malcolm  J. 
Bowman  f and  Yu-Heng  Tseng 


1.  Motivation  and  objectives 

Gulf  Stream  (GS)  separation  near  its  observed  Cape  Hatteras  (CH)  separation  location, 
and  its  ensuing  path  and  dynamics,  is  a challenging  ocean  modeling  problem.  If  a model 
GS  separates  much  farther  north  than  CH,  then  northward  GS  meanders,  which  pinch  off 
warm  core  eddies  (rings) , are  not  possible  or  are  strongly  constrained  by  the  Grand  Banks 
shelfbreak.  Cold  core  rings  pinch  off  the  southward  GS  meanders.  The  rings  are  often 
re-absorbed  by  the  GS.  The  important  warm  core  rings  enhance  heat  exchange  and, 
especially,  affect  the  northern  GS  branch  after  GS  bifurcation  near  the  New  England 
Seamount  Chain.  This  northern  branch  gains  heat  by  contact  with  the  southern  branch 
water  upstream  of  bifurcation,  and  warms  the  Arctic  Ocean  and  northern  seas,  thus 
playing  a major  role  in  ice  dynamics,  thermohaline  circulation  and  possible  global  climate 
warming.  These  rings  transport  heat  northward  between  the  separated  GS  and  shelf 
slope/Deep  Western  Boundary  Current  system  (DWBC).  This  region  has  nearly  level 
time  mean  isopycnals.  The  eddy  heat  transport  convergence/divergence  enhances  the 
shelfbreak  and  GS  front  intensities  and  thus  also  increases  watermass  transformation. 
The  fronts  are  maintained  by  warm  advection  by  the  Florida  Current  and  cool  advection 
by  the  DWBC.  Thus,  the  GS  interaction  with  the  DWBC  through  the  intermediate  eddy 
field  is  climatologically  important. 

Besides  GS  separation  and  mean  path,  realistic  simulation  of  GS  natural  variability 
is  important.  GS  separation  occurs  (by  definition)  when  the  shelfbreak  current  from  the 
Florida  Straits  angles  toward  deeper  water.  If  it  does  not  separate  at  CH,  the  observed 
separation  location,  but  separates  much  later  as  is  common  in  North  Atlantic  Ocean  mod- 
els, it  then  remains  too  close  to  the  shelfbreak  to  allow  northward  (shoreward)  meanders 
to  develop  freely  and  pinch  off  warm  core  eddies  (rings).  The  distance  from  the  mean 
GS  current  to  shelfbreak  must  be  at  least  comparable  to  the  Rossby  radius  of  deforma- 
tion for  this  to  happen.  The  highly  energetic,  nonlinear  turbulent  GS  system  provides  a 
significant  amount  of  geophysical  “noise"  to  the  climate  system.  Such  “noise”  may  affect 
climate  through  fundamentally  nonlinear  stochastic  resonance  (Velez-Belchi  et  al.  2001); 
such  resonance  may  be  associated  with  nonlinear  effects  of  changes  relating  to  freshwater 
sources  and  sinks  (Rossby  and  Nilsson  2003),  changes  of  North  Atlantic  Deep  Water  for- 
mation (Rahmstorf  and  Alley  2002)  and  nonlinear  GS  interaction  with  the  DWBC  seen 
in  results  of  the  present  study.  Such  variability  may  be  natural  in  the  ocean  by  itself. 
For  example,  even  with  annual  cycle  forcing,  constructed  using  observed  annual  cycle 
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surface  wind,  heat  flux  and  freshwater  source/sink  conditions,  there  may  be  significant 
interannual  variability,  as  seen  in  Mediterranean  Sea  modeling  studies  (Fernandez  et  al. 
2003).  Another  example  of  particularly  strong  natural  interannual  variability  is  the  Gulf 
of  Mexico  (Dietrich  et  al.  1997),  which  is  dominated  by  the  Loop  Current  and  its  irregu- 
lar major  eddy  shedding,  independent  of  annual  cycle  forcing,  having  a mean  time-scale 
of  about  270  days. 

Two  recent  major  North  Atlantic  Basin  model  comparison  studies  are:  DYNAMO 
(1997)  and  DAMEE-NAB  (2000).  None  of  the  DYNAMO  models  at  1/3°  horizontal 
resolution  was  able  to  simulate  the  observed  GS  separation  at  CH.  This  failure  using 
only  moderate  horizontal  resolution  was  not  unexpected,  because  of  the  importance  of 
fine-scale  vorticity  dynamics  and  bathymetry  near  CH  (Dengg  et  al.  1996).  Coastal  abut- 
ments such  as  CH  often  lead  to  boundary  current  separation;  examples  are  ubiquitous  in 
nature,  a prime  example  being  the  Black  Sea  whose  coastal  circulation  is  dominated  by 
recirculation  eddies  in  the  wakes  of  coastal  abutments  (Staneva  et  al.  2001).  Besides  iso- 
bath curvature  near  the  CH  abutment,  the  upstream  convergence  of  isobaths  also  plays 
an  important  role  in  GS  separation  (Stern  1998). 

In  the  DAMEE  model  intercomparison  experiment  using  HR  (Hellerman  and  Rosen- 
stein  1983)  climatological  wind  forcing,  all  but  one  of  the  models  failed  to  simulate  GS 
separation  at  CH.  The  exception  was  ROMS  (Haidvogel  et  al.  2000),  which  got  GS  sepa- 
ration at  the  proper  location  without  resolving  the  major  CH  coastal  abutment.  Hurlburt 
and  Hogan  (2000)  report  that  when  sufficient  resolution  is  used  to  address  nonlinear  GS 
separation  dynamics  at  CH  using  the  NLOM  model  (Hurlburt  and  Thompson  1980), 
realistic  GS  separation  and  other  inherently  nonlinear  features  are  well  simulated  using 
HR  winds,  even  though  the  linear  version  of  the  same  model  shows  unrealistic  wind- 
stress-curl-induced  GS  separation  for  11  different  wind  stress  climatologies  including  HR 
winds  (Townsend  et  al.  2000).  The  linear  effects  of  wind  stress  curl  do  not  explain  the 
observed  GS  separation  and  path,  but  instead  give  two  paths,  both  very  different  from 
the  observed  one.  Thus,  the  ROMS  low  resolution  GS  separation  at  CH  is  surprising. 

Herein,  we  adapt  the  DieCAST  ocean  model  to  the  North  Atlantic/Caribbean  Sea/Gulf 
of  Mexico  system  in  order  to  explore  GS  separation  and  ensuing  dynamics.  Base  case 
results  are  discussed,  focusing  on  GS  interaction  with  the  DWBC  during  a 75-year  sim- 
ulation. In  follow  up  to  this  base  case,  closely  related  sensitivity  studies  (Dietrich  et  al. 
2003b)  show  significant  effects  of  horizontal  viscosity  and  diffusivity,  especially  relating 
to  the  thin,  narrow  DWBC  and  its  effects  on  GS  separation. 

The  DieCAST  ocean  model  adaptation  for  the  present  North  Atlantic  study  is  de- 
scribed in  Section  2.  In  Section  3,  the  modeled  North  Atlantic  general  circulation  is 
compared  with  observations,  focusing  on  the  GS  and  DWBC.  Results  before  and  after 
the  establishment  of  the  DWBC  show  dramatic  effects  of  the  DWBC  on  the  GS  separa- 
tion and  path.  Further  discussion  of  the  significance  of  the  DWBC  is  given  in  Section  4. 
Conclusion  and  future  work  are  given  in  Section  5. 


2.  Mlodel  and  Experiment  Design 

We  use  the  z-level  DieCAST  ocean  model  adapted  to  the  North  Atlantic  Ocean/ 
Caribbean  Sea/  Gulf  of  Mexico  system.  The  duo  grid  model  runs  at  about  200  model 
days  per  day  of  computer  time  on  a PC  using  a single  2.0  gigahertz  P4. 

The  model  uses: 

a)  fourth-order  accurate  numerical  approximations  for  all  terms,  except  a conventional 
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Figure  1.  Western  (1/6°  resolution  £ne  grid)  Figure  2.  Eastern  (1/2°  deg  resolution  coarse  grid) 

domain  bathymetry  (depth  is  in  km).  domain  bathymetry  (depth  is  in  km)  showing  shelf- 

break  bathymetry  in  the  shortcircuited  Arctic  Ocean. 


second-order  accurate  hydrostatic  vertical  pressure  gradient,  in  a semi-collocated  control 
volume  framework  (Sanderson  and  Brassington  1998); 

b)  an  incompressibility  algorithm  having  low  numerical  dispersion  associated  with  its 
required  interpolations  (Dietrich  1997) 

The  semi-collocated  grid  avoids  the  large  numerical  dispersion  resulting  from  evaluat- 
ing the  large  Coriolis  term  on  the  conventional  staggered  Arakawa  “C”  grid.  Fourth-order 
accurate  advection  further  reduces  numerical  dispersion.  Thus,  the  model  is  formally  ac- 
curate in  all  significant  respects.  It  is  also  robust  using  very  low  total  (explicit  plus 
numerical)  dissipation. 

The  model  domain  covers  the  North  Atlantic  basin  from  10°  N to  73°  N and  from  97.5° 
W to  0°  W.  To  reduce  the  computation  required,  a duo  grid  approach  is  used;  the  grids 
are  coupled  using  upwind-based  boundary  flux  approximations  as  used  by  Dietrich  and 
Mehra  (1998)  for  the  Santa  Barbara  Channel  nested  in  the  California  Current  system. 
The  western  North  Atlantic,  Gulf  of  Mexico  and  Caribbean  Sea  require  high  resolution 
to  resolve  the  GS  separation  region  and  critical  narrow  straits.  West  of  60°  W,  1/6° 
resolution  is  used;  east  of  60°  W,  1/2°  resolution  is  used.  The  grids  are  fully  two-way 
coupled  each  time  step,  with  one  coarse  grid  cell  overlap  (3x3  fine  grid  cells).  The 
results  are  virtually  seamless  at  the  duo  grid  interface  (Dietrich  2002).  There  are  30 
model  layers,  geometrically  expanding  from  41.6  m thick  at  the  top  to  738  m thick  at 
the  bottom  (maximum  depth  5000  m).  Open  southern  boundary  conditions  are  derived 
from  a one  degree  global  implementation  of  the  DieCAST  ocean  model. 

Figure  1 shows  the  western  domain  bathymetry.  Figure  2 shows  the  eastern  domain 
and  has  lower  resolution  than  Figure  1.  To  avoid  potentially  disastrous  unphysical  vortex 
stretching  in  the  northeast  corner  of  the  eastern  domain,  the  real  bathymetry  is  replaced 
by  a shelfbreak  patterned  after  the  shortcircuited  Arctic  Ocean,  rather  than  using  a 
conventional  vertical  wall  approach. 

To  resolve  the  critical  Caribbean  Sea  passages  and  reduce  computation,  the  duo  grid 
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Figure  3.  Surface  layer  annual  mean  Yashayaev  temperature  (°C)  climatology  in  the  western  domain  GS 
region  (west  of  60  deg  W)  of  the  duo  grid  DieCAST  north  Atlantic  model  (left  panel);  Vertical/latitudinal 
cross-section  of  annual  mean  Yashayaev  potential  temperature  (°C)  climatology  at  70  deg  W,  showing 
highly  nonlinear,  non-diffusive  flattening  of  isopycnals  between  Gulf  Stream  and  Grand  Banks  (right  panel). 


interface  is  placed  just  east  of  the  important  Caribbean  Sea  passages,  and  the  western 
fringe  of  the  Labrador  Sea  (west  of  60°  W)  is  excluded  and  an  idealized  shelfbreak  is 
used  along  60°  W in  the  eastern  domain  (Figure  2),  again  to  minimize  unrealistic  vortex 
stretching.  Surface  boundary  conditions  are  derived  from  monthly  climatology.  The  wind 
forcing  is  by  monthly  average  HR  climatological  winds.  Levitus  climatology  is  used  for 
initial  conditions  and  for  surface  thermohaline  forcing  until  model  year  13,  when  it  is 
replaced  by  an  improved  (less  smoothed)  surface  climatology  (Yashayaev  2002).  The 
annual  mean  surface  temperature  in  the  GS  region  from  Yashayaev  climatology  (Figure 
3)  is  consistent  with  the  observed  GS  separation  near  CH. 

Below-surface-layer  model  temperature  is  restored  toward  climatology  in  a buffer  zone 
20  gridpoints  wide  along  open  northern  lateral  boundaries.  Salinity  is  not  restored,  be- 
cause that  does  not  conserve  salt  material  and  has  no  physical  basis.  Instead,  a 0.2  Sv  (1 
Sverdrup  equals  one  million  cubic  meters  per  second)  inflow  of  freshwater  is  spread  uni- 
formly along  the  open  northern  boundary  to  parameterize  Arctic  Ocean  net  river  inflow 
(Bacon  et  al.  2002)  and  a 0.018  Sv  Gulf  of  St.  Lawrence  freshwater  volume  source  is  also 
specified  (Chapman  and  Beardsley  1989).  Heat  and  freshwater  fluxes  at  the  sea  surface 
are  computed  such  that  the  model  multi-year  mean  annual  cycle  of  surface  tempera- 
ture and  salinity  follow  the  observed  climatological  annual  cycle  (Dietrich  et  al.  2003a). 
This  new  surface  buoyancy  flux  condition  avoids  problems  (reduced  annual  cycle  ampli- 
tudes and  phase  lags,  and  excessive  damping  of  surface  fronts)  attributed  to  conventional 
restoring  (Kill worth  et  al.  2000). 

All  externally  specified  inflows  are  in  the  eastern,  coarse  grid  domain;  these  are  added 
to  the  western  boundary  inflow  from  the  fine  grid  western  domain  and  to  the  integrated 
evaporation-precipitation  (e-p)  to  get  a net  eastern  domain  inflow  (derived  as  described 
by  Dietrich  et  al.  2003a);  then,  the  total  inflow  is  subtracted  uniformly  at  the  southern 
boundary  (about  0.02  Sv  as  it  turns  out)  to  be  consistent  with  incompressibility.  An 
alternate  approach,  not  used  here,  is  to  apply  this  volume  correction  at  the  surface  (e.g., 


Figure  4.  Vertical/latitudinal  cross-sections  at  70  deg  W of  time-averaged  salinity  (left)  and  temperature 

(right)  during  years  0-5. 


in  a global  model  having  closed  lateral  boundaries),  using  the  mean  diffusivity  between 
the  top  two  layers  to  maintain  observed  climatological  annual  cycle  surface  layer  salinity. 
This  approach  is  consistent  with  the  long-term  climatological  state  (near  zero  mean  sea 
level  change)  and  with  the  rigid-lid  approximation  used  in  this  model,  in  which  there  can 
be  no  mean  surface  height  change.  Thus,  the  rigid  lid  is  slightly  porous;  the  divergence 
of  the  barotropic  mode  (vertically  integrated  horizontal  velocity)  is  very  small  but  not 
exactly  zero. 

In  the  western  (1/6°  resolution)  domain,  we  use  constant  horizontal  viscosity  and 
diffusivity  of  10  m2/sec.  Based  on  the  scale  of  the  CH  abutment  curvature  (~  100  km) 
and  GS  velocity  (~  1 m/sec),  this  gives  Reynolds  number  0(10,000)  for  the  primary 
GS  separation  scale.  The  Rossby  number  for  these  scales  is  0(0.1).  Both  domains  use 
near-molecular  background  vertical  viscosity  and  diffusivity,  with  modified  Pacanowski 
and  Philander  Richardson  number  based  mixing  as  described  by  Staneva  et  al.  (2001). 
Thus,  the  inertia  terms  are  strong  in  the  CH  region,  as  necessary  (but  not  sufficient,  as 
a robust  DWBC  is  also  required,  as  will  be  shown  here)  for  nonlinear  inertial  separation. 


3.  Results 

3.1.  Before  and  After  DWBC  Interaction  with  GS  Separation 

The  transient  GS  behavior  during  the  first  five  model  years  is  strikingly  different  from  the 
long-term  behavior  (Figures  4-6).  During  the  first  model  year  (not  shown),  bathymetry- 
induced  strong  inertial  dynamics  lead  to  early  partial  GS  separation  near  CH.  However, 
the  separated  portion  of  the  model  GS  moves  northward  to  the  Grand  Banks  shelfbreak 
during  the  next  few  model  years.  It  then  moves  back  southward  to  its  observed  CH 
separation  point  and  mean  downstream  path  after  about  10  model  years,  and  remains 
close  to  the  mean  observed  path  throughout  years  10  — 75  (end  of  the  present  simulation). 

It  follows  that  model  intercomparisons  of  the  GS  region  before  10  years  of  simulation, 
such  as  done  under  DAMEE,  are  of  limited  value,  especially  when  overly  diffused,  highly 
unstable  (baroclinically)  Levitus  climatology  is  a basis  for  model  comparison,  even  for 
time  averaged  model  results,  as  will  be  seen  below. 

Figures  4 and  5 shows  vertical/latitudinal  cross-sections  at  70°  W of  temperature 
averaged  during  years  1-5  and  years  66  — 75  respectively.  There  is  clearly  a strong 
flattening  of  isopycnal  surfaces  (well  approximated  by  isotherms)  between  the  DWBC  and 
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Figure  5.  Vertical/latitudinal  cross-sections  at  70  deg  W of  time-averaged  salinity  (left)  and  temperature 

(right)  during  years  66-75. 
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Figure  6.  Time  mean  seasurface  height  (cm)  and  700  m depth  velocity  vectors  in  the  1/6°  resolution 
western  domain  during  years  0-5  (left)  and  during  years  66  — 75  (right). 

GS  core  during  the  latter  years  (Figure  5),  as  in  observations  (Figure  3).  Also  during  the 
latter  years,  there  is  a strong  deepwater  GS  jet  and  front  outcropping,  as  in  observations. 

Figure  6 shows  the  time  mean  surface  pressure  with  superposed  700  m depth  time- 
averaged  velocity  vectors  during  years  1 — 5 and  years  66  — 75.  The  color  shaded  contour 
field  is  the  time  mean  surface  height.  The  arrows  are  time  mean  velocity  vectors  at  700m 
depth.  The  mean  Gulf  Stream  front  is  the  boundary  between  the  green  and  light  blue. 
The  slight  northward  bend  into  the  mid  Atlantic  Bight  (a.k.a.  ’’New  York  Bight”)  reflects 
warm  core  eddies  that  separate  from  northern  tips  of  Gulf  Stream  meanders  (mostly  far 
downstream  from  Cape  Hatteras),  some  of  which  slowly  propagate  all  the  way  back  to  the 
Cape  Hatteras  region  where  they  rejoin  the  Gulf  Stream.  Nearly  all  of  the  ’’blue  water” 
in  the  mid  Atlantic  Bight  (a.k.a.  ’’New  York  Bight”)  is  cold  water  that  has  upwelled 
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Figure  7.  Time  mean  pressure  and  velocity  vectors  in  the  1/6°  resolution  western  domain  during  years 
66  — 75  at  levels  1474  m (left)  and  2020  m (right). 


along  the  shelfbreak  between  south  Florida  and  the  mid  Atlantic  Bight  region;  this  water 
recirculates  back  to  the  northern  seas  and  Arctic  Ocean  after  heat  exchange  with  the 
warmer  green  water  that  recirculates  southward  along  with  the  other  ” warmer  colors” 
that  indicate  higher  surface  elevations.  Apparently,  the  DWBC,  indicated  by  the  deep 
velocity  vectors  during  years  66  — 75,  but  missing  during  years  1 — 5,  plays  a big  role  in 
the  observed  mean  GS  path. 

The  DWBC  extends  to  more  than  2000  m depth  in  the  model  results  (Figure  7).  After 
turning  sharply  counterclockwise  near  CH  most  of  the  DWBC  current  travels  eastward 
to  well  off-shore,  before  turning  southward,  crossing  under  the  GS  and  turning  back 
towards  the  coast.  Transient  eddies  are  significant  throughout  the  depth,  and  dominate 
in  the  deeper  levels  where  the  mean  flow  is  smaller. 

These  figures  show  the  dramatic  effect  of  the  DWBC  on  the  GS  and  the  major  flat- 
tening of  isopycnal  surfaces  that  occurs  after  the  DWBC  dense  water  reaches  the  Grand 
Banks  shelfbreak  region.  Dietrich  et  al.  (2003b)  find  that  a realistically  narrower  west- 
ern Labrador  Sea  shelfbreak  significantly  intensifies  the  DWBC,  resulting  in  even  more 
strongly  flattened  isopycnals  than  seen  in  Figure  5,  more  like  Yashayaev  climatology  (Fig- 
ure 3).  Higher  model  resolution  may  further  intensify  the  DWBC,  but  is  not  addressed. 

The  negative  contours  in  Figure  6 intersect  the  U.S.  mainland  coast  nearly  tangentially; 
this  reflects  upwelled  DWBC  water  on  the  west  side  of  the  Florida  Strait,  as  consistent 
with  the  cool  coastal  climatological  temperature  (Figure  3)  and  model  temperature  that 
is  close  to  the  climatology.  This  upwelling  reflects  a secondary  circulation  effect  associated 
with  bottom  friction  and  mixing  of  offshore  GS  jet  momentum  into  the  slower  moving 
coastal  water.  The  strong  GS  jet  effects  dominate  over  local  wind  forcing,  unlike  the 
California  Current  system  (Haney  et  al.  2001). 

Thus,  even  the  nearshore  part  of  the  coastal  current,  which  is  shoreward  of  the  upwelled 
coastal  water  indicated  by  the  negative  contours  in  Figure  6,  separates  almost  completely 
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at  CH  and  forms  the  north  wall  of  the  separated  GS.  As  noted  above,  this  does  not  occur 
as  a linear  response  to  the  HR  winds  used  in  this  study. 

The  explanation  of  the  interesting  strong  GS/DWBC  interaction  indicated  by  Figures 
4-6  is  as  follows.  The  Levitus  climatology  used  to  initialize  the  model  is  unrealistically 
baroclinic  between  the  observed  GS  mean  location  and  the  Grand  Banks  shelfbreak,  re- 
flecting diffusive  smoothing  of  climatological  data.  The  resulting  baroclinic  eddies  trans- 
port potential  density  southward  (mainly  by  northward  heat  transport)  and  flatten  po- 
tential density  contours  in  the  deep  offshore  water,  leading  to  strong  isopycnal  slopes  in 
the  shelfbreak  region  and  much  weaker  slopes  along  the  observed  mean  GS  axis.  Thus, 
there  is  a northward  displacement  of  the  GS  front  to  the  Grand  Banks  shelfbreak  during 
the  first  few  model  years  and  a corresponding  wrong  GS  path.  Only  after  full  establish- 
ment of  the  DWBC  along  the  Grand  Banks  shelfbreak,  which  takes  about  10  model  years 
to  occur,  does  the  front  move  southward  to  its  observed  location.  This  is  the  time  scale 
required  for  densest  DWBC  water  formed  in  the  GIN  Sea  to  reach  the  Grand  Banks. 

Thus,  overly  diffusive  models,  carefully  tuned  to  the  diffuse  Levitus  climatology,  can 
closely  match  this  climatology.  However,  it  is  impossible  to  tune  models  requiring  such 
large  diffusion  to  the  observed  flattened  time  averaged  non-diffused  climatological  isopy- 
cnals  (Figure  3).  This  can  only  be  done  by  a model  such  as  the  present  DieCAST  model 
configured  with  low  diffusion  (Figure  5). 

Thus,  the  DWBC  water  plays  a major  material  role  in  model  GS  separation.  In  a 
limited-area  model  having  specified  DWBC  inflow,  Thompson  and  Schmitz  (1989)  came 
to  a similar  conclusion.  The  role  is  analogous  to  the  inner  wall  of  rotating  annulus  exper- 
iments, except  the  cold  water  source  is  by  advection  of  the  dense  DWBC  water  rather 
than  by  conduction  at  the  wall.  The  GS,  in  turn,  plays  a role  analogous  to  the  outer 
wall  of  the  rotating  annulus  experiments.  Amplitude  vacillation  in  those  experiments  is 
accompanied  by  flattening  of  interior  isotherms  due  to  baroclinic  instability  (Pfeffer  et 
al.  1974),  which  by  its  very  nature  results  in  such  flattening.  This  also  characterizes  the 
atmospheric  index  cycle. 

3.2.  Other  Significant  Results 

A large  northern  recirculation  gyre  occurs  between  the  GS  and  the  Grand  Banks  shelf- 
break at  depth  > 700  m.  This  gyre  entrains  and  recirculates  some  of  the  DWBC  water. 
Figure  8a  shows  its  far  north  source,  the  East  Greenland  Current,  flowing  over  the  Den- 
mark Strait,  which  preconditions  even  deeper  water  formation  in  the  Labrador  Sea.  On 
the  east  side  of  the  Denmark  Strait,  a fragment  of  the  North  Atlantic  Gyre  flows  into  the 
GIN  Sea  as  observed.  Notably,  such  two-way  flow  separated  by  a strong  front  occurs  in 
spite  of  the  coarse  1 /2°  resolution  used  in  the  eastern  domain  of  the  present  model.  The 
sill  overflow  quickly  becomes  a quasi-balanced  bottom  density  current,  as  the  Coriolis 
terms  quickly  arrest  any  initial  downslope  flow  (away  from  shore)  to  give  a quasi-balanced 
along-isobath  flow  with  deep  water  and  lower  pressure  to  the  left  (looking  downstream); 
a slow,  frictionally  induced  cross-isobath  (away  from  shore)  flow  toward  deeper  water 
occurs  in  the  bottom  boundary  layer. 

Figure  9 shows  the  time  mean  rms  velocity  deviation  from  the  time  mean  velocity 
at  700m  depth  in  the  Denmark  Strait  region,  where  dense  sill  level  water  spills  “over 
the  dam  (sill)”  into  a bottom  density  current.  This  shows  that  the  deep  eddy  activity 
is  concentrated  in  the  DWBC  region  (along  the  shelfbreak),  where  significant  vortex 
stretching  is  driven  by  energy-releasing  sinking  dense  DWBC  water  along  the  shelfbreak 
and  by  locally  strong  interaction  with  GS  water  near  the  northeast  coast  of  the  Flemish 
Cap  (45  deg  W,  47  deg  N).  The  vortex  stretching,  as  the  potential  energy  is  released, 
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Figure  8.  (a)  Vertical/longitudinal  cross-section  of  time  mean  latitudinal  velocity  across  Denmark  Strait 
in  the  coarse  1/2°  resolution  model  eastern  domain  (left  panel);  (b)  Time  mean  surface  height  and  velocity 
vectors  at  700  m depth  in  the  Denmark  Strait  region  (right  panel). 


results  in  turbulent  mixing  of  the  East  Greenland  Current  (Arctic  Ocean  source)  water 
with  the  North  Atlantic  Gyre  water.  The  saddle  point  structure  of  the  time  average 
surface  height  (Figure  8b)  reflects  the  sill  bathymetry,  and  also  shows  the  mean  watermass 
paths  leading  to  turbulent  interaction.  Thus,  Figures  8-9  indicate  inertial,  turbulent 
interaction  between  the  East  Greenland  Current  and  North  Atlantic  Gyre  south  of  the 
Denmark  Strait,  even  with  the  relatively  coarse  1/2°  resolution  in  this  region.  The  sill 
overflow  is  an  important  part  of  the  thermohaline  circulation,  and  the  associated  locally 
strong  mixing  is  important  to  watermass  dynamics. 

While  some  of  the  DWBC  upwells  on  the  north  side  of  the  GS  separation  point,  some 
of  it  continues  southward.  Some  of  the  southward  branch  upwells  on  along  the  shelfbreak 
between  CH  and  the  Bahamas.  The  cool  coastal  temperatures  in  both  the  Yashayaev 
climatology  and  the  model  results  in  this  region  reflect  this  upwelling.  Horizontal  advec- 
tion  by  the  Florida  Current  from  the  south  that  warms  the  coastal  region  is  balanced 
by  vertical  advection  and  mixing.  Such  vertical  circulation  is  forced  by  a combination 
of  bottom  drag  and  entrainment  of  coastal  water  into  the  strong  GS  current.  This  indi- 
rect circulation  intensifies  the  GS  front,  thereby  affecting  development  of  eddies  through 
baroclinic  instability,  especially  after  GS  separation.  Thus,  in  this  narrow  coastal  region, 
the  upwelled  undercurrent  water  is  entrained  into  the  Florida  Current  from  the  south; 
this  is  a region  of  significant  watermass  mixing.  Unlike  the  wind  forced  coastal  jet  of  the 
California  Current  system  (Haney  et  al.  2001),  local  wind  effects  may  be  secondary  in 
such  an  energetic  current  system. 

These  mean  results,  calculated  on  a modern  PC,  are  more  realistic  than  the  recent 
state-of-the-art  results  by  Chao  et  al.  1996  despite  requiring  less  overall  resolution.  Using 
realistic  HR  winds  which,  unlike  the  winds  they  used,  do  not  force  separation  at  CH 
(Townsend  et  al.  2000),  the  mean  path  of  the  model  GS  is  very  close  to  the  observed  mean 
path,  starting  with  a realistic  trajectory  from  the  observed  separation  point.  Further 
downstream,  the  separated  upwelled  coastal  water  bifurcates  from  the  GS  core,  similar 
to  the  two-branch  separated  Gulf  Stream  in  the  linear  response  to  HR  winds  noted  by 
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Figure  9.  Time  mean  rms  velocity  deviation  from  the  time  mean  velocity  at  700  m depth  in  the 

Denmark  Strait  region. 

Townsend  et  al.  (2000),  but  is  also  influenced  by  interaction  with  the  New  England 

Seamount  Chain,  as  indicated  by  the  animation  from  Figure  10  (shown  in 

ftp : //ssc . ere  .msstate . edu/pub/mehra/nab_duo/OMO/year69-70 .AVI . gz) . 

Figure  10  shows  a sequence  of  the  surface  currents  during  model  year  70.  A major 
cold  core  eddy  (often  called  a ring)  pinches  off  the  southern  tip  of  a southward  GS 
meander  around  day  85.  Around  day  100,  a warm  core  pinches  off  a northward  meander. 
Later,  another  warm  core  develops  as  a result  of  GS  interaction  with  the  New  England 
Seamount  Chain.  There  is  a deep  northward  secondary  flow  generated  by  blocking  of  the 
deep  GS  flow,  which  also  steers  the  northern  GS  surface  current  fringe  northward  into 
a developing  warm  core  ring,  as  discussed  by  Hurlburt  and  Thompson  (1984).  This  GS 
bifurcation  is  a prevailing  feature  in  model  results,  as  indicated  by  Figure  6.  The  GS 
straightens  out  to  its  observed  mean  path  after  these  two  major  eddy  pinchoffs.  These 
results  are  consistent  with  observations  of  transient  small-scale  features  in  and  around 
the  GS.  The  results  of  Chao  et  al.  (1996)  do  not  show  such  fine-scale  features. 

Several  times  per  year,  buoyant  warm-core  eddies  pinch  off  northward  meanders  of  the 
GS  and  drift  westward  and  southwestward  back  through  the  Middle  Atlantic  Bight,  as 
also  seen  in  Figure  10.  They  become  imbedded  in  the  Middle  Atlantic  Bight  shelfslope 
current  and  ultimately  entrained  back  into  the  GS  near  CH  (e.g.  Bowman  and  Duedall 
1975).  However,  the  volume  of  slope  water  entrained,  including  its  eddies,  is  only  a 
fraction  of  the  total  0(100  Sv)  GS  transport.  Full  separation  of  the  warm  GS  water 
occurs  at  CH,  including  a small  volume  of  much  cooler  DWBC  water  upwelled  along  the 
shelfbreak  between  the  Florida  Strait  and  the  Middle  Atlantic  Bight. 

Also  in  the  Figure  10  sequence,  the  Gulf  of  Mexico  Loop  Current  penetrates  northward 
and  pinches  off  a major  warm  core  eddy.  The  Dry  Tortuga  cyclone  plays  a significant 
role,  as  do  other  cyclonic  frontal  eddies  that  originate  along  the  Yucatan  shelfbreak, 
some  of  which  merge  into  the  Dry  Tortuga  eddy  when  the  Loop  Current  is  well  extended 
into  the  Gulf  of  Mexico.  Some  frontal  cyclones  also  pass  through  the  Florida  Strait. 
Before  the  major  warm  core  eddy  pinchoff,  the  Loop  Current  takes  the  shape  of  a square 
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Figure  11.  Root  mean  square  deviation  of  surface  height  from  years  66  — 75  time  mean  in  western  1/6° 

resolution  domain. 

with  rounded  corners  due  to  frontal  eddy  effects.  The  Loop  Current  takes  a right  angle 
turn  into  the  Florida  Strait  after  eddy  shedding.  Westward  propagating  anticyclones 
prevail  in  the  Caribbean  Sea,  and  there  is  a semi-permanent  cyclonic  recirculation  in  the 
southwestern  Caribbean  Sea  corner.  These  details  are  consistent  with  observations  and 
earlier  modeling  results  (e.g.  Dietrich  et  al.  1997). 

Figure  11  shows  rms  seasurface  height  anomaly  derived  from  the  model  results.  In 
general,  the  rms  value  is  30  — 45  cm  along  the  separated  GS  axis  all  the  way  to  the  New 
England  Seamount  region  near  the  eastern  boundary  of  the  1/6°  resolution  grid.  Thus, 
vigorous  transient  eddies  populate  the  modeled  GS  region  with  root-mean-square  surface 
height  anomaly  amplitudes  comparable  to  those  in  nature  (DYNAMO  1997,  Figure  8.6). 
The  rms  anomaly  in  the  Gulf  of  Mexico  region  is  also  close  to  the  observed,  as  in  earlier 
Gulf  of  Mexico  studies  with  DieCAST  (Dietrich  1997;  Dietrich  et  al.  1997).  The  rms 
value  in  the  coarse  resolution  domain  east  of  60°  W (not  shown)  is  smaller,  but  there  is 
a local  maximum  of  about  20  cm  along  the  shelfbreak  at  the  south  edge  of  the  Flemish 
Cap.  Coupling  to  low  resolution  coarse  grid  domain  may  cause  reduced  rms  values  near 
the  eastern  boundary  in  the  fine  grid  domain. 
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4.  DWBC  Interaction  with  Methane  Hydrates:  Potential  for  Global  Warming 

Accurate  DWBC  simulation  is  important,  because  it  not  only  affects  the  GS  path  as 
shown  in  Section  3,  but  the  DWBC  path  also  affects  potential  catastrophic  release  of  vast 
amounts  of  methane  gas  from  methane  hydrates  stored  on  the  shelf  slope  below  depths 
where  methane  freezes  (about  4°  C,  depending  also  on  pressure-depth).  It  is  thought  that 
this  methane,  after  escaping  to  the  atmosphere,  may  oxidize  and  add  huge  amounts  of 
CO2  (Katz  et  al.  1999),  a greenhouse  gas  that,  unlike  the  water  vapor  by-product,  does 
not  precipitate  out,  thus  leading  to  global  warming. 

As  suggested  by  Lai  (2003),  strongly  exothermic  biogeochemical  processes  occuring 
within  the  methane  bubbles  generated  in  the  ocean  depths  (Kennett  et  al.  2000)  may 
significantly  warm  the  DWBC,  lead  to  positive  feedback,  including  by  DWBC  recircula- 
tion in  the  northern  seas  and  Arctic  ocean,  and  affect  thermohaline  circulation  including 
Arctic  and  northern  seas  stratification  through  annual  cycle  ice  melting/freezing.  These 
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processes  are  coupled  to  the  GS.  The  nonlinear  GS/DWBC  system  dynamics  may  lead  to 
accelerating  methane  release  and  sudden  global  warming.  As  further  noted  by  Lai  (2003), 
this  possibility  demands  careful  study,  due  to  observed  recent  global  warming  (IPCC 
2001)  and  oxygen  depletion  (Joos  et  al.  2003)  in  the  ocean  depths  near  the  methane 
hydrate  levels.  A major  concern  is  that  anthropogenic  perturbation  of  the  system  may 
lead  to  even  stronger  catastrophic  warming  than  the  about  10°  C Arctic  temperature 
rise  that  occurred  in  10  — 20  years  at  the  end  of  the  last  glacial  maximum  (Rossby  and 
Nilsson  2003). 

The  GS  path  affects  the  location  and  intensity  of  North  Atlantic  deep  water  formation, 
which  occurs  mainly  in  the  northern  seas  and  Arctic  Ocean  and  feeds  the  DWBC  after 
spilling  over  the  GIN  Sea  sills.  Its  changes  would  alter  the  location  and  temperature  of  the 
DWBC,  possibly  augmenting  the  potential  catastrophic  greenhouse  gas  induced  warming 
noted  above.  Modeling  the  DWBC  water  is  a big  challenge  relating  to  GS  separation. 
Dense  fluid  in  a thin,  narrow,  bottom  density  current  must  travel  from  its  high  latitude 
source  region  to  the  Grand  Banks  shelfbreak  without  excessive  numerical  dissipation 
and  dispersion,  in  order  to  have  sufficient  intensity  to  realistically  affect  the  energetic  GS 
separation  and  ensuing  dynamics;  this  travel  takes  0(10)  years. 

The  dense  narrow  thin  DWBC  is  not  resolved  by  the  climatological  data  used  to  initial- 
ize the  model.  Thus,  it  must  develop  spontaneously  in  the  model.  This  dense  shelfslope 
water  forms  in  the  northern  seas,  and  takes  more  than  a year  to  reach  the  Grand  Banks 
from  the  Labrador  Sea,  and  0(10)  years  to  reach  there  from  the  GIN  Sea.  In  order  for 
it  to  arrive  at  CH  with  sufficient  intensity  to  properly  affect  GS  separation  and  path,  it 
must  not  be  overly  diluted  by  numerical  dispersion  and  dissipation. 

5.  Conclusion  and  future  work 

Results  presented  herein  show  that  one  may  simulate  GS  dynamics  and  separation  from 
the  coast  realistically  in  a marginally  eddy-resolving  simulation  by  using  a model  having 
a)  low  total  (physical  plus  numerical)  dissipation;  b)  fourth-order  accurate  advection 
with  its  associated  low  numerical  dispersion;  c)  further  numerical  dispersion  reduction 
in  the  interpolations  between  grids;  and  d)  fourth-order  accurate  horizontal  pressure 
gradient.  Nonlinear  dynamics  are  involved  in  GS  separation  and  it  thus  seems  necessary 
to  resolve  the  major  CH  abutment.  However,  although  bathymetry-induced  strongly 
inertial  dynamics  may  favor  separation  near  CH,  our  results  indicate  that  is  not  sufficient; 
a robust  DWBC  is  also  required.  A robust  DWBC  requires  that  the  DWBC  water  not 
be  dispersed  by  numerical  errors  during  the  0(10)  years  to  travel  from  its  source  region 
to  the  Grand  Banks  shelfbreak  region.  The  importance  of  the  DWBC  is  enhanced  by  its 
interaction  with  methane  hydrates,  which  may  lead  to  global  warming. 

It  would  be  interesting  to  nest  a finer  grid  in  the  GS  separation  region  to  shed  more 
light  on  the  detailed  dynamics  of  the  GS  separation,  focusing  on  the  vertical  circulation 
in  the  region  of  GS  separation  and  how  this  circulation  connects  the  GS  to  the  DWBC. 
The  significance  of  simulating  a realistic  GS/DWBC  system  lies  partly  in  its  potential 
impact  on  climate  through  fundamentally  nonlinear  stochastic  resonance  effects  that  are 
known  to  occur  in  many  nonlinear  systems  (Rahmstorf  and  Alley  2002). 
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